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Abstract 

We continue to study an infinite-parametric family of gauge theories with an arbitrary function 
of the self-dual part of the field strength as the Lagrangian. The arising one-loop divergences are 
computed using the background field method. We show that they can all be absorbed by a local 
redefinition of the gauge field, as well as multiplicative renormalisations of the couplings. Thus, 
this family of theories is one-loop renormalisable. The infinite set of /3-functions for the couplings 
is compactly stored in a renormalisation group flow for a single function of the curvature. The flow 
is obtained explicitly. 


1 Introduction 

The purpose of this paper is to compute the one-loop beta-function for an infinite-parametric family 
of gauge theories introduced in [Tj. The Lagrangian of these theories is given by an arbitrary function 
of the self-dual part of the field strength. In paper [1] we have performed a perturbative study of these 
theories, and analysed the arising scattering amplitudes. We have also computed the beta-function 
for the first non-trivial coupling describing deformations away from Yang-Mills. It was noted that this 
family of theories can be expected to be renormalisable at the one-loop level, by a local redefinition of 
the gauge field, as well as by multiplicative renormalisation of the coupling constants. In this paper 
we use the background field method to compute the arising divergences, and establish the one-loop 
renormalisability. We also extract the full one-loop beta-function for an infinite number of coupling 
constants parameterising our theories. 

Given that the calculation that leads to the /3-function is quite technical, it is worth stating the 
result already in the Introduction. The family of theories we study is described by the following 
Lagrangian 

£ = M 4 /(FWM 2 ). (1) 

Here M is a mass scale introduced so as to make all couplings dimensionless. The quantity Ff IN is 
the self-dual part of the field strength for the gauge field, a is the Lie algebra index, and M, N are 
unprimed spinor indices, see below for our spinor index conventions. The function / is an arbitrary 
gauge- and Lorentz-invariant function of Ff IN . In case this function is analytic near zero, it can be 
expanded in powers of the argument. The first non-trivial term is then the YM Lagrangian, and 
the next term, cubic in the field strength, describes the first non-trivial chiral ’’deformation” of YM 
theory. It is clear that a theory of type (JT]) is specified by an infinite number of coupling constants, 
coefficients of the power series expansion of the function /. Apart from the YM term quadratic in the 
curvature, all other interactions described by (H|) are power-counting non-renormalisable and contain 
inverse powers of M 2 . 
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We compute one-loop divergences using the background field method. We work in Riemannian 
signature. The one-loop renormalisability of © means that, after a local field redefinition, the re¬ 
maining one-loop divergences are taken care of by counter terms of the type already contained in (fTj). 
Thus, the one-loop running of the whole infinite set of the coupling constants can be encoded as the 
running of the dimensionless function / 


df(x) 
d log fi 




( 2 ) 


where /3f is some gauge- and Lorentz-invariant function of the dimensionless self-dual part of the 
curvature x a AB := F AB /M 2 . The result of our calculation is 


Pf( x ) = 
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Here ( f) aAB and (/") 


tf\abABC D 


are the matrices of the first and second derivatives of the function / 
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and (f") 1 is the matrix inverse to f". We also use the notation 

^ab . riaeb / rf\ab s~iaeb( £f\e 

x AB ■— u x ABi \J JAB ■— ° \J ) ABi 


where C abc are the Lie algebra structure constants. 

Note that the result ([3]) is a homogeneity degree zero function in /, as it should be. As a further 
check, we note that for Yang-Mills 
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/ y m(x) = j-^ab) , 


(5) 


and thus 


(/ym )ab — 2^2 x abi ((/ ym ) 1 )abcd — 2 g 2 5 ab e A (c e \B\D)- 


A simple computation then gives 


p™ 


11 c 2 

6(47t) 2 


Mb) 2 , 
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where C 2 is the quadratic Casimir. This gives the correct running of l/4g 2 . 

When both the Lagrangian (|T]) and the /3-function © are expanded in powers of x AB , one can 
read off an infinite set of beta-functions for the couplings stored in /. In the main text we carry out 
this exercise for the first non-trivial coupling parameterising the T 3 interaction. We reproduce the 
result obtained in [4J by a different method. Some further comments on the interpretation of © are 
contained in the Discussion section. 

Let us give an outline of how m is computed. All details are given in the main text. The 
straightforward application of the background field method to © runs into a difficulty. The problem 
is that the second-order operator that arises by linearising m is not of Laplace-type, even after 
gauge-fixing. We alleviate the problem by passing to the first-order formulation, by ’’integrating in” 
an auxiliary field. The linearisation of the resulting Lagrangian then gives a first-order operator that 
turns out to be of Dirac-type. Its square is of Laplace-type, which makes the well-developed heat-kernel 
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technology for operators of such type applicable. This allows for a straightforward, even thought at 
times technical, computation of the divergent parts of the relevant regularised determinant. 

The organisation of this paper is as follows. In the next Section we reduce the problem to a 
calculation of a certain determinant, which is then analysed in Section [3j We apply the result obtained 
to the problem at hand in Section [H This is where the result (J3]) is obtained. An application to a 
2-parameter family of theories is considered in Section [5] where the result obtained in [T] is reproduced. 
We conclude with a discussion. 

2 Setup of the calculation 

2.1 Non-Laplace type operator 

Before we start the calculation, let us explain why we need to do more than just apply the standard 
background field method to (0). The problem with the straightforward application of the method is 
that the linearisation of the Lagrangian m around an arbitrary background 

C ( 2 ) = l{fy bA ^ D {da) a AB (da) b CD + M\fT AB [[a , a]] a AB (7) 

gives a non-Laplace type operator. Here a a AB is the gauge field perturbation, and the Lagrangian is 
expanded to second order. The matrices f, f" are the already encountered 0 matrices of the first 
and second derivatives of the function /. The other quantities are defined as follows 

(da) AB := 2d( AA ra B j A , [[ a ;°]]AB := C abc a b AA ta c B A . (8) 

Here d is the covariant derivative with respect to the background gauge field. It is clear that the 
second order differential operator appearing in (0) is not of the Laplace type because the indices of the 
derivatives d AA > contract either with a spinor index of the gauge held perturbation or with a spinor 
index of the matrix of second derivatives. Thus, we have a second-order operator A acting on Lie 
algebra-valued one-forms a“ with the leading second order part of the form 

(Aa)“ = M ab ^d u d p a b a + ..., (9) 

for a certain matrix M abpi,pa . This matrix is not proportional to rf' p , and thus the operator does not 
reduce to a Laplace-type one with the leading part M abpa d 2 a b a for some matrix M abp,a . Thus, the 
well-developed heat-kernel technology for Laplace-type operators is not applicable in our case, at least 
not directly. 

2.2 The first-order formulation 

To apply the heat-kernel technology to 0 we would need to develop the associated expansion for 
non-Laplace type operators, which is possible, but not easy. We will choose another route. Thus, we 
will first reformulate our theories in the first-order form by ’’integrating in” an auxiliary field. The first 
order operator arising in this context turns out to be of Dirac-type, with square given by an operator 
of Laplace-type. This will reduce the problem at hand to the well-studied Laplace-type operators. 
Let us consider the following action principle 

S[A,B] = J B aMN Flj N - M 4 V(B/M 2 ), (10) 

where F^ N is the self-dual part of the curvature of A a MM ,, and B^ N is an auxiliary field. The field 
equations that follow from (jlOli are 

d A E B\ E = 0, F a MN /M 2 = ( V'Tmn , ( 11 ) 
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where 


is the matrix of first derivatives of the potential function V(B) 


( V' 


l\aMN 


dV 

d(B aMN /M 2 )' 


(12) 


We can solve the second field equation for B/M 2 in terms of F/M 2 , and substitute the result to (1101) 
to obtain a theory of type (flj). This identifies V(B/M 2 ) as the Legendre transform of the function 
f(F/M 2 ). So, we have an equivalent formulation, at least classically. One can also check explicitly 
that the exact equivalence continues to hold at the one-loop level, see below. 

The usual Yang-Mills theory can be described in this language, and corresponds to the choice 


M a V ym (B/M 2 ) = g\B a MN ) 2 . 


(13) 


Indeed, in this case, integrating the auxiliary field B out we get B = (1/2 g 2 )F and 

5 YM [A] = ±J (FS, n ) 2 , (14) 

which is the usual Euclidean signature YM action. 

In general, held equations (jlll) imply a second-order differential equation for the connection. First, 
one solves the second equation for B^ N as a function of F^ N . This is an algebraic task. In the second 
step one substitutes the solution into the first equation, and obtains a second order PDE for A a MM ,. 
It is easy to see that this is the same PDE that is obtained directly from the Lagrangian (JT|) . 


2.3 Linearisation 

We now want to compute the one-loop /1-function for (jlOl) using the background held method. As we 
shall see, the associated problem can be reduced to a Laplace-type operator. We first linearise the 
action (1101) around an arbitrary on-shell configuration of helds B,A. We denote the fluctuations by 
b, a. The resulting quadratic in the fluctuations Lagrangian is 

£ ( 2 ) = 2 b aNM d MM >a a N M ' + C abc B aMN a b MM ,a c N M ' - \{V")t NPQ b aMN b bPQ . (15) 

Here cLaa' is the covariant derivative with respect to the background connection, and the matrix 
cn mnpq i s that of second derivatives of the potential function V(B). The reason why in the first 
term we wrote b aNM with the order NM of the spinor indices will become clear below. 

Before we proceed with our analysis, let us check that the functional determinant arising from (1151) 
is the same one that arises from (|7|) . Thus, let us integrate out the perturbation of the auxiliary field. 
The field equation for b a MN reads 

{da) a MN = (y") < M N PQb bP Q ■ (16) 

Given that the function V is the Legendre transform of / we have 

(( y") 1 )mnpq = U"TmnpQi (17) 

and thus the solution of (fTHl) is 

Kin = U") a MNPQ{do) bPQ ■ (18) 

Using the background field equation B^ N = (,/Omao we then see that (fl5l) with b a MN integrated out 
is precisely (p7|) . This establishes the one-loop equivalence of the first- and second-order formulations. 
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2.4 Gauge-fixing 

The kinetic term in csi) is degenerate, as it is built from the operator that sends the linearised 
connection field a a AA , with its 4 x diiri(g) components into an object n) M w hh 3 x dim(g) 

components. We thus need to gauge fix in order to be able to invert the kinetic term. 

The gauge-fixing fermion can be taken to be the usual one in the context of Yang-Mills theories 

T = 2 c a d%,a a M M ' - 2M ab c a h b , (19) 


where c a is an anti-ghost, h a c is an auxiliary field, and we have allowed ourselves to introduce some 
matrix M ab instead of the gauge-fixing parameter. This can, in principle, depend on the background 


field B a MN . 


The BRST variation of (fTUj) is 


sT = 2(?d%,d M M 'c a + 2 h a c df v a a M M ' - 2M ab h a c h b . 


( 20 ) 


The first term here gives the kinetic term for the ghosts, while the second term should be added to 
the bosonic part of the Lagrangian, and removes the degeneracy discussed above. 

A particularly convenient way to combine the last two terms in (1201) with m is to introduce a 
new component to the b a MN field 

~ h aNM = b aNM + M ^ ( 21 ) 

The new field b aNM is no longer MN symmetric, and the two indices become distinct. This explains 
why we wrote b aNM in the kinetic term in () 15 D . With the new field introduced, the first term in (1151) 
plus the second term in (I20j) become simply 2b aNM dMM ,a % M ■ Further, he last term in (1201) combines 
with the last term in (1151) by introducing 

C V")mnpq = ( V")mnpq ~ M ab e M NepQ. (22) 

The full gauge-fixed Lagrangian then becomes 

£ ( 2 ) = 2 b aNM d M M'a a N M ' + f abc B aMN a b MM ,a c N M ' - \(V")t NPQ b aMN b bPQ , (23) 

which is the same as m , but now for the field b. From now on we will omit the tilde from the 6-field 
for brevity. 


2.5 The operator arising 

The operator in (1231) is of the type 


D = 


B 


(24) 


where 


Q : V®S. -> F® 5+, <f : V ® S + -> V ® S. 
are the usual chiral Dirac operators and 

AeEnd(Y®5+), BGEnd(F®5_). 

Here S± are the spaces of unprimed and primed spinors. In the case of (1231) the bundle V is given by 

V = g ® S+ , (25) 


5 


and it is the first index N of b aNA1 that is considered as taking values in the S+ part of (1251) . Further, 
in our case the endomorphism A is simply 

Aab = A^anB = ~2 (y")MANB- ( 26 ) 

Finally, the endomorphism B is 

Ba’B' = B^ IN eA'B' = C aeb B e MN €A'B' ■ (27) 

This can be written in more abstract notations as 

End(V®S-) 3B = B-1 S _, B € End(H), (28) 

with 1 s_ being the identity operator on S-. In our case 

We now proceed with the calculation of the (divergent part of) the determinant of (|24l) . and return 
to the problem at hand afterwards. 


3 Calculation of the determinant 

3.1 Problem setup 

Consider the Dirac-type operator (1241) . with A,B being some symmetric endomorphisms, so that the 
operator is self-adjoint. We assume that B is of the form (1281) . This operator acts on columns 

beV®S + , a£V®S-, 

where V is some vector space. In other words, the operator (1241) acts on E-valued Dirac spinors, and 
it is assumed that the chiral Dirac operators <ft, (jf know how to act on E-valued sections. 

We would like to compute the divergent parts of the regularised determinant of the operator (1241) . 
One way to do this computation would be to eliminate either a or b fields, and consider the resulting 
second-order operator. We already know that if we integrate out b we get a non-Laplace type operator. 
It is easy to check that what one obtains is exactly the operator in dTj). Thus, eliminating b is not a 
good idea. Another possibility is to eliminate a, and obtain a second order differential operator for b. 
It is easy to see that this has a better chance to work, because the arising operator is of Laplace type. 
However, in the process of integrating out a one needs to invert the matrix B, which is possible, but 
cumbersome. 

Instead of trying to reduce the problem at hand to a second-order Laplace type operator by 
integrating one of the two fields a, b, let us consider the square 

2= ( #* + A 2 A? + ?B \ 

V B(ft* + fA + B 2 J 

It is clear that the arising second order operator is of Laplace type, and hence the divergences of 
interest are easily computed using the heat-kernel techniques. 

To perform the computation we will need to rewrite D 2 as a multiple of 

- VV,, + E, (31) 

where T> is some covariant derivative and E is some endomorphism. We will compute both T> and E 
in the following subsections. 


(30) 
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3.2 Spinor index conventions and Lichnerowicz formuli 

To do calculations with (1241) we will use spinor index notations. Thus, our conventions for the actions 
of fy, $* are as follows. Let G S + ,X A i G S-. Then 

OX)A = d AA ,\ A \ (, pX)A' = d A ' A \ A , 


where d AA ' is the usual covariant derivative operator in spinor notations. 

Let us use the spinor index notations to compute the squares (j)(f and We expect these to be 
multiples of the corresponding Laplacians plus curvature terms, so let us work out what appears. We 
have 

(<ft<ft*\) A = d AA 'd A E Xe = —~d M M'dM M Aa + - F a e Xe , (32) 

{(j>*(j)X') A i = d A / E d EE 'X E = —~d M M 'dM M X A r + -F a , e ,X e . 


Here F AE ,F A / E i are the self and anti-self-dual parts of the curvature of defined by the following 
formula 

‘■^djid,/ X — F^ u A, 

and the same for the action on X'. Then the following decomposition of T) Jt , G End(H) is valid 


Fmnvnn' = -Fmn^M'N' + - Fm'N'Cnm, 


(33) 


so that 


M' 


Fmn — F MM i N 
Equalities (1321) show that 2 D 2 is an operator of type (1311) . 


£ M 'N'=r M'MN' ■ 


3.3 Computation of the connection 

To compute the connection that makes 2D 2 into an operator (1311) we first rewrite the off-diagonal 
elements in a suggestive form 

2= ( -d»d li + F + 2A* -2C% + 2(?5) A 

V -2C*% + 2(Q*A) -d^d^ + F' + 2B 2 ) ' 1 j 

Here 

F G End(H (g> S + ), F' G End(V r (8) 5_) 
are the self- and anti-self-dual parts of the curvature viewed as operators, and 

C* : V <g> S+ -A V <g> S_ 

are operators acting as 

{c„ X') A = c A A ,x A \ (c*x) a > = c; a a ,x a . 

The components of C^, C* are given by 

C M m ,A a> = £ A ’M’{A + B1 S+ ) AM i C* M m’ A A' = £m'A'(A + Bt s+ ) MA . ( 35 ) 

Finally, the derivatives of B,A appearing in (1341) are as follows 

($B) : V®S- -» y®5+, (?*A) : F <g> S+ -> 
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with the explicit expressions being 


am \') A = (d AA 'B) x ai 

This immediately tells us that 


OP A)\) A . = (d A ' A A A E )\ E . 


-D — I dp 
M 1 c; d. 


and 


E = 


F + 2A 2 + 2 {m + ^C'v) 


2 (<f A) + (d»C*) F' + 2 B 2 + C;C 
After some simplifications the endomorphism becomes 

F - 2 (AB + BA) - 2 B 2 • l s (0(5 - A)) 


E = 


\A-B)) 


F' - (Tr(A 2 ) + Tr(A)5 + B Tr(A)) • l s _ 


( 36 ) 


(37) 


(38) 


where we introduced new maps 

- A)A')a = (d MA '(B ■ 1 5+ - A) a m )X a ', (f(A - B)\) A , = (« d A , A (A - B ■ 1 s+ )a E )\e- 


3.4 Curvature 

The curvature of is given by 


T _ 9T) r, _ ( F^ + C.Ct-CvC; 2d [lJ C„ ] 

- iu^u v] - ^ 2 d [M c* F, w + c;c v - c;c „ 

3.5 The heat-kernel formula 

The quantity we want to compute is, see e.g. [2] 

° 4 = ( 4 ^ 1 + 


(39) 


(40) 


3.6 Off-diagonal terms 

First derivatives of the matrices A, B enter only via the off-diagonal terms. These off-diagonal terms 
do not mix with the diagonal when the trace is taken. It is thus pertinent to analyse these clearly 
separable derivative terms first. 

Consider the contribution from the term 6 E 2 first. Squaring the off-diagonal term produces two 
diagonal terms that become the same under the trace in V. Thus, we can write for this contribution 

12TV y ( d M A ,{A - B • 1 S+ )AMd NA ' {A - B • 1 S+ )tv A ) • (41) 

The other contribution is from F 2 . Again, there are two equal terms and we get 

8T \ v (d [tl C v]AA ,d^C*^ AA ') . (42) 

While (1411) and (1421) are both quadratic in first derivatives of A, B, the space-time indices are contracted 
differently. We would like to bring the expression (|42j) into the form similar to (1411) , plus contributions 
without derivatives. This is achieved by a sequence of identities. 



The idea is first to rewrite the contraction of two anti-symmetric tensors in (|42l) as the contraction 
of the corresponding self-dual parts, plus a term containing The necessary identity is 

2 d [p C u] d^C^ = 2 (d M , {M C N) M ')(d^C* N ^ N ') - e^d p C u d p C*, (43) 

where we suppressed the spinor indices of C P ,C* for brevity. We can integrate by parts in the last 
term, and modulo a total derivative write 


e^d,c u d p c; = \t ,wpa CpF, w C* = \c™Fmn ° C* NM ' - ]-C MM ' F w N , o Cf' 


2 - — ~ c t 2 ~ ivi ' ~ Lv± 1 v 2 ~ - tvi iv m ■ (44) 

Note that here the curvature acts on the quantity C* that takes values in End(E), not in V. To get 
this identity we have used the following spinor expression for the 4-dimensional e-tensor 

e livpa _ e M'N' e R' S' e MR e NS _ ^M'R' ^N'S' ^MN ^RS 


as well as (1331) . In our case the last term in (PHI) does not contribute because both C p ,C* are 
proportional to the primed spinor metric, which contracting with the symmetric Fm’N' gives zero. 
Thus, reintroducing the spinor indices on C AA , C* AA we get for (1421) 

- 8 Try (d A , [M X AN) d A '( M X N ^ + 4 Try ( X A M F MN ° X NA ) , (46) 

where we denoted 


Xab — {A + B • 1s + )ab, (47) 

and in the last term the curvature acts on an object in End(E &S+), which is denoted by the o symbol. 

The final identity that we need is one reducing (1411) to an expression of the type appearing in (1461) . 
This is again achieved by an integration by parts. We have 

d%X AM d A ' N X N A = -X AM d%d A ' N X N A = -X AM d A ' N d%X N A - X AM F mn o X n a (48) 

= d A 'NX A Md A ' M X NA - X AM F mn o X n a . 

But then we have 

Aa'nXam = d A '(MXAN) ~ d%X ae■ (49) 

Substituting this into (1451) we get 

d%X AM d A ' N X N A = d A ' ( m X AN )d A> ( M X N ^ A - ld%X AM d A ' N X N A - X AM F mn o X n a (50) 

and thus 

d A > {M X A N)d A ' {M X N ) A = | d%X AM d A ' N X N A + X AM F mn o X n a , (51) 

which is the identity we need. 

Collecting all the above, we see that the total contribution of the off-diagonal terms is given by 
12 Try (d M A ' {A - B • 1 S+ )AMd NA ' (A - B ■ 1 S+ )n A ) (52) 

-12Try [d M a’ {A + B ■ ls+) A Md NA (A + B ■ ls + )jv A ) 

-4 Try (Xa M F MN °X na ), 

where we kept the matrix (1471) in the last term. Thus, in the terms containing the derivatives of A, B 
only the cross-term remains and we get for the contribution from the off-diagonal terms 

24 Try (d AA 'Tr(A)d AA ' Bj - 4Try (X A M F MN X NA ) + 4Try (X MA F MN X A N ) , (53) 

where Tr(^4) = Ae e ■ We see that this is the only part of A whose derivative appears in the result. 
Note that we also expressed the action of the curvature on an object in End(C (g> S+) as a sum of two 

terms. In each of the last two terms in (1531) both F. X are viewed as objects in End(C <8> S+). 
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3.7 Diagonal terms 


There is no derivatives of A, B coming from the diagonal terms. Working their contributions is simply 
a tedious algebra. We start with the E 2 contributions. The upper diagonal term contributes 

Try ( -F MN F mn + 4 (A MN B + BA M n)F mn + 8Ti-(A 2 )B 2 + 8 A M N BA N M B + 16Tr (A)B 3 + 8.(54) 

The lower-diagonal term gives 

Try (- F min ,F m ' n ' + 2Tt(A 2 )Tx(A 2 ) + 4Tr(A 2 )(Tr(^4)S + BTr(A)) + 4Tr(A)BTr(A)B + 4(Tr(A)) 2 B 2 ^J .(55) 


To work out the F 2 contributions we first write down what one obtains for the maps on the 
upper-left and lower-right corners. We have 


Fmn(-A B - 2A' yl ( M X A r)' B 


f-M'N’ + -^Fm'N'£NM£A 


B 


for the first diagonal element of the matrix of curvatures and 


-^Fmn + X AN ) 


tM'N'tA' + 


1 


-F M 'N' e A' B — Fr(X 2 )e A i^ M ie N i) 


B' 


£NM 


(56) 


(57) 


for the second. 

Squaring and taking the trace of the first gives 

Try (f MN F mn + F m , n ,F m ' n ' - 4 X MA F MN X A N + A(X AM X n b + X AN X M B )X B M X NjC ) . (58) 
With a bit of algebra the last term here can be rewritten as follows 

Try (4(X AM X n b + XanX m B )X b M X na ) = Try (4Tr(X 4 ) - 8Tr(X 2 )Tr(X 2 )) . (59) 

Thus, overall, the square of the upper-left diagonal element gives 

Try (f MN F mn + F m , n ,F m ' n ' - AX ma F MN X a N + 4Tr(X 4 ) - 8Tr(X 2 )Tr(A 2 )). (60) 

The trace of the square of the lower-right diagonal element gives 

Try (F mn F mn + 4X a M F MN X na + 2(X M A X AN + X N A X AM )X MB X B N ) (61) 

+Try (F M ' N ,F M ' N ' - 6Tr(X 2 )TV(X 2 )) , 

where the second line is from the square of the anti-self-dual M'N' symmetric part. The last term in 
the first line can be brought into a more canonical form. Overall, after some algebra we get 

Try [f MN F mn + F M ' N 'F m ' n ' + 4 X A M F MN X NA - 4Tr(X 4 ) - 4Tr(A 2 )Tr(X 2 )) . (62) 


The overall contribution from F 2 is obtained by adding (1601) and (l62lh We get 
Try (2F MN F mn + 2F M i N iF m ' n ' - 4 X MA F MN X A N + 4 X A M F MN X NA - 12Tr(A 2 )Tr(X 2 )) . (63) 
We can now substitute the matrix X. We have 

Try (Tr(X 2 )Tr(X 2 )) = Try (^Tr(A 2 )Tr(A 2 ) + 2Tr(A 2 )(Tr(A)B + BTv(A)) (64) 

+2Tr(A)BTr(A)B + 2(Tt(A)) 2 B 2 + 4T y(A 2 )B 2 + 8Tr (A)B 3 + 4 B'^ . 

Finally, adding (1341) and (1331) . multiplying the result by 6, and adding to (1331) we get 

Try ( - 4 F MN F mn - 4F M 'N'F m ' n ' + 24 (A MN B + BA MN )F MN (65) 

-4 X ma F MN X a n + 4 X A M F MN X NA + 48 A m N BA n m b) , 
with all other terms nicely cancelling away. 
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3.8 The final result 

We can now combine the contribution (1531) from off-diagonal terms with (1651) to get the final result 
for the relevant heat kernel coefficient 

af = ^ ^ J Try (24 d AA fTr(A) d AA ' B - 4 F MN F MN - 4 F M , N ,F M ' N ' (66) 

+24 {A MN B + BA MN )F mn + 48 A m N BA n m b) . 


4 Calculation of the /^-function 

With the hear-kernel coefficient (1661) at hand we can return to our problem of computing the /3-function 
for (HI). We now use (1661) with (f26l) . (1291) . 


4.1 On-shell conditions 

The integrand in (|66l) is not of the form (HOD of the action we started from. In particular, the first 
term in (1661) contains derivatives of the auxiliary field B. It is easy to see that this term goes away for 
the case of the usual YM theory, when the matrix A A b is a constant. However, for a general theory 
of the type we consider this derivative term is non-zero. 

To eliminate the term in (I66D that contains unwanted derivatives, and also to be able to identify the 
/3-function, we need to go on-shell and assume that the background satisfies its field equations. So, let 
us derive some consequences of the background held equations (fTTT) . We already mentioned that the 
held equations can be viewed as giving rise to second order PDE’s on the connection. However, we can 
also obtain some second order equations for the auxiliary held. Indeed, applying another derivative 
to the first equation in (fill) we get 

0 = d A A 'd A , E B a BE = - l -d MM 'd MM 'B% A - l -C aeb F% E B b BE . (67) 

We note that the second term on the right-hand-side is automatically AB symmetric as a consequence 
of the second equation in (HID . Indeed, the AB anti-symmetric part is a multiple of 

C aeb F e MN B bMN = C aeb {V') e MN B bMN = 0, (68) 


where the hrst equality is using the second of the held equations, and the last equality is a simple 
consequence of the gauge-invariance of V(B). 

It is now clear that using the on-shell condition (|67j) we can replace the derivative term in ()66D by 
a term that contains no derivatives of B. To this end, it is convenient to write (I67D as follows 


0 = d A A 'd A , E B ab E = --d MM ' 


dMM'B BA - -C 


aec 77>e E rich 
b A &be 


^ s-ibec 77 e E fyac 
— —O Jr a d 


BE • 


(69) 


Now, introducing the matrix F^J N € End(H) 


F] 


ab 

MN 


riaeb t? 
.= G r ; 


MNi 


(70) 


we can rewrite (I69D in the matrix form as 


— d MM dMM'BBA = F a e B E b — BebF a e . 


(71) 
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4.2 On-shell simplifications 

Let us do some rewriting of the first term in (RUl) . Substituting (17TT) we get 

- 24Tr(7l) (AS) (. F a e B eb + B b e F ea ). (72) 

Note that we have added the symmetrisation on Tr(j4) AB , which is for free because the quantity in 
the brackets is automatically AB symmetric (as a consequence of field equations). 

Let us now consider the first term in the second line in (1661) . It is given by 


^{Aam B nB b a + B A R Abm a n)F 


B 


-iMN 


(73) 


where all Lie algebra indices are suppressed and simple cyclic index contraction is assumed. We now 
assume that no gauge-fixing has been done on the A-matrix and it is symmetric in its pairs of indices. 
There is no harm in this assumption because the gauge-fixing part of A anyway drops from the final 
result, as can be explicitly checked. The final trick is to do some exchanges of the spinor indices. For 
the first term this gives 


1 


A A m b nB b a = A nm b aBb a + e A NA E M Bh B b a = -Fmn + Ti(A)m b B B n, 

where we have used the background field equation in the form 

1 


(74) 


B mb A%, cd = ~I& D . (75) 

We have now spelled out the Lie algebra indices for clarity. For the second term we get similarly 

(76) 

Thus, overall, we can rewrite (|73l) as follows 


B a A B m n = B a A BM n + £m B a A be n = -Fmn + B m Ay(A) bn . 


2 \F AB F AB + 24Tr(A) A " (B b *F ea + F b *B ea ), 


\AB 


(77) 


where we have relabelled some indices. It is now clear that the first term (1721) cancels some of the 
term just analysed. Indeed, we can further rewrite (1771) as 

24 F ab F ab + 24TV(7l)( AB ) ( B b e F ea + F b e B ea ) + 24Tr (A) E E B MN F MN , (78) 


where we have used 


f)MN 77« _ t?MN f> 

& -TMN ~ £ &MN, 


(79) 


which is the already discussed consequence of the field equations. Now the second term in (178|> cancels 
with (1721) . Thus, overall, on-shell the relevant heat-kernel coefficient becomes 

1 1 


af = 


8 F AB F AB + 2AA ab ab B mn F MN + 48A am BN B b c A C n DM B d a 


(4tt) 2 24 

where we divided the previous result by 2 to obtain the coefficient for the original operator D. 
We can now add the ghost contribution 


(80) 


a 


ghost _ 


l l 

( 4^)2 12 


f ab f 


AB 


to get 

D o ghost 

a 4 - 2al = 


1 1 


(47t) 2 6 


FarF ar + 6 A ab ab B mn F MN + V2A am rn B b u A cn um B n 


AB tdMN 


BN j 


C 


DM 


(81) 


(82) 


which the final /3-function quoted in Q. In the formula (j3j) we have also wrote A, B explicitly in 
terms of the matrices of first and second derivatives of the function /. We note that this formula is 
only valid for the original matrix A that is symmetric in its both index pairs. In other words, this is 
already a gauge-invariant formula from which the gauge-fixing parameters dropped out. 
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5 Application to a 2-parameter family 

5.1 2-parameter family of theories 

Let us consider a gauge theory containing just a single additional interaction term as compared to the 
standard YM 


£ - 4^2 ( f mn ) 2 + 


a 


3lg 2 M 2 


z^iabc rpaB rpb C rpc A 
° r A ^ B £C • 


(83) 


The sign in front of the first term is as appropriate for the Euclidean action. When the coupling 
constant is absorbed into the connection one obtains ag in front of the new term, which is a convenient 
parameterisation of the new interaction, as we shall see. It is our goal to obtain the /3-function for ag. 

To be able to use the result of the computation above, we first need to put (1831) into the first order 
form. This is achieved by the following Lagrangian 


C = B 


aMN tp a 


F a MN - 9 {Bmn) + 


a(2g 2 ) 2 

3 M 2 


C abc B a A B B b B c B c c a + 0{B 4 


Indeed, varying with respect to B we obtain 

a{2g 2 ^ 2 


f mn ~ ^9 2 B% n + 


M 2 


C abc B bB B c EN + 0{B A ). 


We can find B AIN from here perturbatively, as a series expansion in F. We get 


1 

27 2 


r>a _ rpCL 

a MN ~ 7T2 B MN 


a 


C abc FlfF c EN + 0(F 6 ). 


2 g 2 M 2 

Using the Jacobi identity we can rewrite this in the following form convenient for the latter 


r)db _ ^ Tpdb 

&MN — ~ ~ B i\/f j\j — 


a 


pae E rp€.b 


2 g 2 1 MN g 2 M 2 { m N)E- 

Substituting this into (1841) we get (1831) . plus higher order corrections. 

5.2 The matrix of second derivatives 

The matrix of second derivatives of the potential function is then as follows 


1 


2 ag 4 


AmaNB ~ TmANB ~ ~9 F 1 ^M(N^AB) °{ e N(MB A ) B + e B(MB C A ) N )- 

Using spinor identities we can manipulate the last term into a more convenient form 

AmANB = ~2 9~3 ab (eMNtAB + £MB£An) - J^P^MnB'ab + e A bBmn )> 

where we have used the notation (1291) . The full trace is now easy to compute 

= -3 g 2 5 ab . 


Aab MA 
Fl a 


(84) 


(85) 


( 86 ) 


(87) 


l MA 


( 88 ) 


(89) 


(90) 
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5.3 Computation 


Let us now substitute all the above into the final result (1521) . Using C apq C bpq = C 2 d ab the first term 
becomes —C 2 {F AIN ) 2 . For the second term we get 


6(_% ( 2^ F ™ n ~ 


a 


g 2 M 


rpae E rpeb 
2 r {M r N)E 


pbaMN = 9 C 2 {F ^ n ) 2 + 


9aC 2 

M 2 


(jabc jpa B jpb C 


C > 


(91) 


where we used 


(japb (jbqc (jcra 



(92) 


The computation of the last term in (182(1 is a bit harder, but the result is 1/3 times (1911) . Adding 
everything up we get 


a 


D 

4 


— 2 a 


ghost 

4 


c 2 1 

( 47 t ) 2 6 


11 ( f mn ) 2 + 


F^OL /-<a6c T?aB rpb C 77c A 
JpF t A t B t c . 


(93) 


5.4 Interpretation 


The result (1931) directly gives the /3-functions for the coefficients appearing in the Lagrangian (1831) . 
Indeed, we can immediate write 


d /J_\ = C 2 11 
d log /i \4<7 2 / (47 t) 2 6 


d ( a \ = C 2 
dlogg \3\g 2 M 2 J (47r) 2 M 2 


(94) 


From here we get the familiar 

dg 11 g 3 K 

dlogg 3 (47r) 2 ’ 


(95) 


which is unchanged in our deformed family of theories. The second equation gives, together with the 
/3-function for the YM coupling 


d(ag) ag 3 K 
dlogg (47t) 2 ’ 

which coincides with the result obtained in [I] by a different method. 


(96) 


6 Discussion 

In PQ we have motivated out interest in the family of theories dTJ) by a possibility that this family 
of theories may be renormalisable in the sense of effective field theory, i.e. after all possible held 
redefinitions are taken into account. The main result of this paper is that this is indeed the case at 
the one-loop level. We have also explicitly computed the arising renormalisation group flow ([3]). 

Perhaps a good analogy indicating why our result ([3]) is interesting is to the non-linear cr-model in 
2 + e dimensions. This is a model with the Lagrangian 

m=j 9iMix))d tl 4> l {x)d IJ ,(j^{x), (97) 

where (p l {x) is a map from a two-dimensional surface to some target space with metric gij. Because 
in 2 dimensions the held (p l is dimensionless, all powers of the held are allowed as renormalisable 
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interactions, and so the above field theory is renormalisable. It is a famous result [3| that the /?- 
function for the infinite number of coupling constants stored in the metric (jij ((f)) is given by the Ricci 
tensor Rij , so that the renormalisation group flow is the Ricci flow. 

In this paper we have computed the one-loop renormalisation group flow ([3]) for theory (jT|) . The 
result is a flow in the space of scalar functions of the self-dual part of the field strength, and is in a 
sense an analog of the Ricci flow for (1971) . It would be very interesting to understand a geometrical 
interpretation of ([3]) , if any, but we leave this to future work. 

The most important open question is whether the renormalisability continues to hold at higher 
loops, or is lost at some sufficiently high loop order. Indeed, any held theory is renormalisable once 
all terms compatible with the symmetries are included into the Lagrangian. This is renormalisability 
in the sense of effective held theory. However, the effective held theory Lagrangians with their inhnite 
number of terms of ever increasing number of derivatives are too complicated to explicitly compute the 
associated renormalisation group how. We are guided by a vision that there may exist four-dimensional 
renormalisable effective held theory models that come with an inhnite number of couplings, but that 
are still sufficiently simple so that explicit calculations of the renormalisation group how in the infinite- 
dimensional space of couplings are possible. The class of theories studied in this paper is an attempt 
to construct such a model. 

Our result that the theories (JT|) are one-loop renormalisable says little about what is to happen 
at higher loops. Indeed, there are examples of power-counting non-renormalisable theories that are 
renormalisable at one loop, but this property is lost at higher loops. One such example is General 
Relativity. Indeed, it is a famous result [4] that GR is one-loop finite with zero A, and one-loop 
renormalisable with A/0. At two loops this property is lost [5]. In the case of GR there is a single 
coupling, and so perhaps this is not a good analogy with the infinite-parametric models considered 
here. Another, more appropriate example is that of scalar field in 4 dimensions with an arbitrary 
potential V((f>). It is an easy exercise to check that this theory is one-loop renormalisable, with the 
one-loop /3-function for V(<f>) being a multiple of ( V") 2 . However, this property is lost at two loops. 

In the case of the theories considered here, our justification for the renormalisability hopes is the 
” collapsing” property of perturbation theory explained in details in [T]. This property implies that 
there is at most a single derivative on internal Feynman diagram lines, and leads to improved UV 
behaviour. In the considered in this paper first derivative formulation this single derivative property 
is realised automatically, because in this formulation the vertices contain no derivatives, and the 
propagator is the inverse of a first-order operator, which thus can be written as a single derivative of 
an inverse of a second-order operator. Thus, it is again the ”collapsing” property of our model that 
allowed us to perform the one-loop computation and obtain an analog of the Ricci flow for (1971) . but 
now in the context of a power-counting non-renormalisable theory. We hope that this property will 
also make the higher loop behaviour nice (renormalisable), but at present we have no calculations to 
support these hopes. 

The question of whether the renormalisability continues to hold at higher loops is related to 
the question whether and in what regime one can trust the flow ([3]). Certainly, even if a model is 
renormalisable, the one-loop flow can only be trusted if higher loop corrections are negligible, i.e. in 
the regime of validity of the perturbation theory. So, if our model continues to be renormalisable at 
all loops, we would just need to compute (or estimate) the associated /3-function corrections to see 
in what regime the one-loop flow can be trusted. However, in case renormalisability is lost at some 
loop level, one needs to add new operators to the tree level Lagrangian. The presence of these would 
affect the one-loop flow, as well as the very ability to compute such a flow. So, at present we do not 
understand how to establish the regime of validity (if any) of a one-loop flow for a theory that becomes 
non-renormalisable at some higher loop order. And this is why hope, possibly over-optimistically, that 
there are models of the type considered here that exhibit all-loop renormalisability and allow for an 
explicit computation of the arising renormalisation group flow. 
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If the flow ([3|) can be trusted, at least in some regions in the space of couplings, the most interesting 
question is where this flow takes one in the ulta-violet. The first term in Q is non-positive, while 
the last term is non-negative. Thus, there can in principle exist functions / which give the vanishing 
/3-function. These can be IR or UV fixed points. Since the theory at such a fixed point must be scale 
free, the corresponding function / must be a homogeneous degree 2 function of its argument, so that 
the dependence on M drops out. However, given that the number of arguments of the function / is 
in general as large as 2dim(G) — 3, there is a vast range of possibilities even if the requirement of 
homogeneity is imposed. At present we are unable to say anything about the existence of fixed points 
other than the one corresponding to free theory that is realised in the usual YM theory. 

If there are non-trivial fixed points, the next question is about the dimension of the critical surfaces. 
For a given fixed point this is the surface of all the trajectories that lead to it. If there are fixed points 
with a finite-dimensional critical surface, one has a realisation of the asymptotic safety scenario of 
Weinberg, see [6] for a recent discussion. What is most interesting about our result ([3]) is that it now 
allows for a very concrete study of these questions, for the first time in the four-dimensional setup. 

It is pertinent to compare our approach to that followed in the modern asymptotic safety literature, 
see e.g. [7]. In this approach, a version j8j of Polchinski’s exact renormalisation group equation [9] is 
used to compute and study the flow. To compute anything in practice, a truncation in the space of 
effective actions must be taken, and then one studies this truncated flow. In the case of our model, 
because of renormalisability, no truncation is necessary, at least in the presently studied one-loop 
setup. Another important difference between this work and the approach reviewed in [7] is that 
in the latter divergences other than logarithmic are taken seriously, and one studies the flow for 
dimensionful couplings. In contrast, in this work we used the dimensional regularisation, kept only 
the logarithmic divergences and obtained the flow for dimensionless couplings. The /3-function we 
computed is independent of the gauge-fixings used at intermediate stages (this has been used as a 
check of the final result). What we obtained is a flow in the space of observable coupling constants, 
with the latter being operationally defined from e.g. on-shell gluon scattering amplitudes. 

Let us further remark that our result (j3j) may be interesting even if the hopes about all-loop 
renormalisability are too optimistic. Another possible scenario could be that there are some smaller, 
perhaps finite dimensional, subspaces in the space of theories of type dH) such that if one is inside the 
subspace the renormalisation does not take one out. Indeed, the usual YM theory is precisely such a 
one-dimensional subspace. It is possible that one can find other subspaces of this type by studying the 
flow Q. This would be interesting, and would give new gauge theories in four space-time dimensions. 

The results of this paper also give support to the conjecture stated in ms that a certain infinite- 
parametric class of four-dimensional gravity theories is (one-loop) renormalisable. In fact, the struc¬ 
tures that arise in the context of these gravity theories are quite analogous to what we have encountered 
here. Similarly to what we have seen for the case of gauge theories CO), a straightforward application 
of the background field method for gravity theories m leads to a non-Laplace type operator, see El- 
One lesson from the calculation performed in this paper is that an efficient strategy for dealing with 
such operators is via the first order formalism. Work on applying these ideas to gravity is currently 
in progress. 
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